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Magnetic Resonance Imaging (MRI) can cliaracterize and discriminate among tissues using their 
diverse physical and biochemical properties. Unfortunately, submicrometer screening of biological 
specimens is presently not possible, mainly due to lack of detection sensitivity. Here we analyze the 
use of a nitrogen- vacancy center in diamond as a magnetic sensor for nanoscale nuclear spin imaging 
and spectroscopy. We examine the ability of such a sensor to probe the fluctuations of the 'classical' 
dipolar field due to a large number of neighboring nuclear spins in a densely protonated sample. We 
identify detection protocols that appropriately take into account the quantum character of the sensor 
and find a signal-to-noise ratio compatible with realistic experimental parameters. Through various 
example calculations we illustrate different kinds of image contrast. In particular, we show how to 
exploit the comparatively long nuclear spin correlation times to reconstruct a local, high-resolution 
sample spectrum. 

I. INTRODUCTION 

Physical tools have historically facilitated advances in biology; notable examples are X-rays crystallography, DNA 
sequencing, microarrays techniques, and, above all, microscopy in its various forms. Extending Nuclear Magnetic 
Resonance (NMR) to the micro- and nano-scale promises to become another leading resource in the microscopist's 
toolbox: Unlike any other technique, NMR is unique in allowing the generation of images with different information 
content. Multidimensional high-resolution spectroscopy is today routinely used in the liquid and solid state to unveil 
complex molecular structures, and this capability could prove groundbreaking if samples having sub-microscopic 
dimensions could be efficiently probed. Unfortunately, these features cannot be fully exploited at present because 
NMR lacks the sensitivity essential to high-resolution screening. The origin of this limitation is twofold: First, in 
"conventional" NMR the signal to noise ratio (SNR) is proportional to the nuclear magnetic polarization of the 
sample, which represents only a small fraction of the attainable maximum 10^'* for protons in a 14 T magnet at 
300 K). Second, Faraday induction is a poor detection method since, even with maximum polarization, the minimum 
number of spins needed to induce a measurable signal is comparatively large. 

Although experiments performed at lower temperatures and/or higher fields can partly mitigate these problems, 
other more efficient detection techniques have recently been proposed. One strategy is to use the spin associated to 
a single nitrogen- vacancy (NV) center in diamond as a local magnetic field probe. [Tl [5]. The operating principles of 
this approach closely mimic those of an atomic vapor magnetometer [3], where the applied magnetic field is inferred 
from the shift in the Larmor precession frequency. Owing to the exceptionally long coherence times of NV centers — 
exceeding 1 ms at room temperature in ultra-pure bulk samples [3] — detection of 3 nT over a measurement time of 
only 100 s has been experimentally demonstrated. [5^ Further, a NV center within a diamond nanocrystal attached to 
an AFM tip was recently used to image a magnetic nanostructure with 20 nm resolution. |6j 

Here we focus on applications of a NV center mounted on a scanning probe for monitoring adjacent nuclear 
spins in an external, infinitely-extended organic sample. Rather than detecting single nuclear spins — an extremely 
challenging goal — we focus on the case where the NV center interacts with large ensembles of nuclear spins localized 
over effective volumes of ^(10-50nm)'^. This regime lends itself to a simplified description that simultaneously takes 
into consideration the quantum nature of the sensor — the NV center — while relying on a classical description of the 
long-range dipolar fields induced by the nuclear spin ensemble. Similar to prior magnetic resonance force microscopy 
experiments [7], our strategy exploits the small dimensions of the effective sample to probe the 'nuclear spin noise' — 
i.e., the statistical fluctuations of the nuclear magnetization — rather than the magnetization itself. An important 
consequence is that, unlike traditional MRI, spatial resolution is not due to strong magnetic field gradients but is 
rather determined by the distance between the NV center and the sample. Assuming a very small external magnetic 
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field we determine the conditions required for 2D nuclear spin imaging at (or near) room temperature, and show 
them to be compatible with realistic experimental parameters. Further, we show that, in addition to determining 
the local nuclear spin density, this strategy allows one to explore different kinds of contrast mechanisms (nearly a 
requisite when imaging, for example, densely protonated organic/biological systems). In particular, we show how to 
reconstruct the local nuclear spin correlation function and, from it, a spatially-resolved nuclear spin spectrum. 

The paper is organized as follows. First, we briefly review the operating principles of NV-center-based magne- 
tometry, more explicitly identify the effective size of the sample being probed, and lay out our detection protocol. 
Subsequently, we describe different modalities of 'nuclear spin noise' detection and determine in each case the limit 
signal-to-noise ratio. Finally, we discuss image contrast and localized nuclear spin spectroscopy and conclude with 
some model calculations. 

II. SPIN-NOISE MAGNETOMETRY WITH A SINGLE NV CENTER 

The negatively charged nitrogen-vacancy center in diamond is an impurity comprising a total of six electrons, two 
of which are unpaired and form a triplet ground state with a zero-field splitting Dgs— 2.87 GHz. In our calculations 
we assume the presence of a small magnetic field Baz (^ 10 mT) collinear with the crystal field (which, in turn, is 
oriented either along the [111] axis or its crystallographic equivalents). Though non-mandatory, the auxiliary field 
lifts the degeneracy between the \ms = ±1) states, thus allowing one to selectively address only one of the two possible 
transitions, e.g., between jm^ = 0) and jm^ — 1). 

When a green laser (532 nm) illuminates the NV center, the system is excited into an optically active triplet 
state; subsequent intersystem crossing produces a dark, singlet state that preferentially relaxes into \ms = 0). Almost 
complete optical pumping of the ground state takes place after a ^ 1/iS illumination, thus allowing us to model the 
initial density matrix of the NV center -for practical purposes, a two-level system- as 

P(0) = |0)(0| = i(/ + a.), (1) 

where / denotes the identity operator and cr^ is the Pauli matrix. Because intersystem crossing is allowed only if 
excitation takes place from |ms = 1), the fluorescence intensity correlates with the population of the spin state. We 
model the 'measurement' operator as 

M = a|0)(0|-|-6|1)(1| = ^{a + b)I+^{a-b)a,. (2) 

In Eq. a and b are two independent, stochastic variables associated with the total number of photons collected 
during the measurement interval (~ 300 ns) and characterized by Poisson distributions qa{k) = a^e~°'/kl and qb{k) — 
P'^e~^ /k\, with k integer. Due to the branching ratio into the dark singlet level, the averages over several measurements 
a = (a) and /3 = (b) are substantially different (a = 1.5/3) and thus provide the contrast necessary to discriminate 
the sensor spin state. 

Fig. la schematically shows the basics of our detection protocol: spin initialization and a selective tt/2 microwave 
pulse are followed by a period At of free evolution in the presence of an unknown, nuclear-spin-induced magnetic 
field Bnz. Preceding optical readout, a second tt/2 pulse, shifted by a phase 9 relative to the first pulse, partially 
converts spin coherence into population differences. In the rotating frame resonant with the chosen transition, the 
density matrix describing the NV center is given by 

p{At) = i(/-(7,sin(0 + 0)+a,cos(0 + ^)), (3) 

where (j) = JJ^* JcBn (t) dt denotes the total accumulated phase due to the nuclear field and 7e is the electronic 
gyromagnetic ratio. As in any other magnetometer-based strategy, the goal of a measurement is to extract the value 
of (j) and, from it, valuable information on the magnetic field. 

Before considering the constraints deriving from the quantum character of the sensor, we describe the magnetic field 
generated by the nuclear spin ensemble. In an experimental setup where the NV center scans an infinitely-extended 
sample film, the electronic sensor spin and the nuclear spins are coupled via long-range dipolar interactions. Given 
that in the rotating frame resonant with the sensor spin only components of the nuclear field parallel to the z— axis 
need be taken into consideration, we find 

BftT =Bj,z = Y, {/ (r-) m« + g (r,) (tt^^^ • f,) z} . (4) 
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FIG. 1: (a) Basic diamond- based magnetometry pulse sequence, (b) With the NV center at the reference frame origin, the 
grayscale indicates the relative contribution to field fiuctuations from spins in a uniformly dense film, (c) In units of the relative 
radial coordinate s/d, the upper set of curves shows a cross section of the graph in (b) (black curve) and the corresponding 
integral (white curve). The grey curve shows the effective spin noise 'density' &jv (s, z) (see text). For comparison, the lower set 
shows the same curves but for the average field at the NV center. Note that the integral (dashed white curve) decays slowly 
to zero as a result of negative contributions from spins far from the center. 



Here f {r) — {^io/4:Trr^) (Scos^ i? — l) and g{r) — (3/xo/47rr^) cosi^ are functions of the distance = |ri| of the z-th 
nuclear spin to the NV center and di is the angle between the position vector and the z— axis; /io is the magnetic 
permeability of vacuum, and mi^ {rn"^ ) denotes the component of the corresponding nuclear magneton m^*^ parallel 
(perpendicular) to the 2;-axis. We will consider the situation where the distance d between the sensor and the surface 
is of order ~ 10 nm or greater. We also assume that nuclear spins are dense (i.e., no nuclear spin can be singled 
out). In this regime, the NV center interacts with a large number of protons -exceeding 10^ in most organic samples- 
and thus exerts a negligible back-action on the sample system. Each nuclear spin can be described classically via a 
stochastic, ergodic variable featuring first and second moments (m) and (m^), respectively. 

To see that detection of the time- dependent fluctuations of the nuclear field — rather than the field itself — better 
suits our purpose, let us consider the case of a uniformly magnetized film and assume, for simplicity, that the normal 
to the sample surface coincides with the z-axis. Using Q we write the time-averaged field acting on the sensor as 



{Bm) I {/ (f) (m,) + g (r) {{m^) ■ f)} dV, (5) 

•'p J Film 

where we have transformed the sums into volume integrals via the correspondence — >■ / ^ with Vp representing 
the volume of the 'primitive cell' associated to a single nuclear spin. From symmetry considerations, we observe 
that the second term in ([5| cancels out. This is also the case for the first term -in agreement with the classical 
magnetostatics result outside a thin, infinitely-extended, uniformly-polarized film- but here a more subtle balance 
between contributions from spins close and far away from the sensor is responsible [H]. The latter is shown in Fig. Ic 
where we plot / (f) (and its integral) as a function of the (normalized) radial coordinate s on the sample plane; within 
each thin slice of thickness dz, long-range, weaker contributions from more numerous spins far from the sensor exactly 
cancel the field created by spins contained within a central disk (of diameter comparable to the sensor-slice distance) . 

The concept of spin noise detection capitalizes on the spontaneous fluctuations of the nuclear spin magnetization 
in a small volume. To more quantitatively identify the sample volume within the film, consider the special case of 
a uniformly distributed, infinitely-extended sample and calculate the nuclear field variance Ai3^. Starting from Eq. 
Q and in the limit of ([s]) we find 

ABl^{Bl)^^ [ {{f{r)f{ml) + {g{r)sindf{ml)]dV, (6) 

J Film ^ ' 
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where we assumed (m^) = ('ti^) = Using cylindrical coordinates for convenience, we plot in Fig. Ic the 'spin 

noise density', bN{s,z)dsdz = ^ Iq^ ^ (/(s, z))^ {ttiD + {g{s, z)sj\J + s^) (m^) S(iiy9| dsdz. While spins far 

from the sensor have a non-negligible contribution, fluctuations of the nuclear field at the NV center are dominated by 
spins approximately contained within half a sphere of radius comparable to the sensor-surface distance d. Comparing 
with the prior results, we conclude that fluctuations selectively highlight spins close to the sensor -as opposed to 
'distant' spins- not because the resulting average field is stronger but because, being less numerous, the relative field 
variance is larger. 

A practical upper limit for the NV center-sample distance d stems from the fact that the amplitude of the field 
fluctuations decreases sharply with the sensor-sample distance: Assuming a sample with spin density pN ~ ^/Vp, we 
find 

ABn - CpomNpT/d''^^ (7) 

with C a constant of order ^1/20 obtained from integration of Eq. and m^vthe nuclear magneton. For example, 
in the case of an organic system with proton density pjv ~ 5 x l(F°m^"^ and assuming d ^ 200 nm, we obtain 
ABn ^ 2.5 nT, a value approaching the sensitivity limit of a room-temperature, diamond-based magnetometer ^[S]. 

We note that detection of the average magnetization within the 'active' volume -as opposed to magnetization 
fluctuations -is conceivable if the contribution to the total field from spins outside this volume has been canceled. [5] 
In this case the nuclear field Bat at the NV center site has the approximate value 

Bn ^ DpomNp^P, (8) 

independent of the sensor-surface distance. Here P — mpfBA/{2kBT) is the nuclear Boltzmann polarization at 
temperature T and D is a, constant of value ~ 1/6. Comparing Eq. ([?]) and ([s]) we find the criterion for spin noise 

dominance, d < {kBT/ (rriNBA))'^^^ Pj^^^^ ■ For example, if we take as a reference the case in which the protonated 
sample (pn '-^ 5 x lO^^m^'^) has been polarized to the equivalent of a magnetic field Ba—W T at room temperature, 
we have d <500 nm. 



III. SENSITIVITY LIMITS 



Having identified the source and magnitude of the field fluctuations at the sensor site, we now turn our attention to 
the general problem of using a quantum object — the NV center — to gather information on the fluctuating ensemble 
of sample spins. The average fluorescence in the presence of the nuclear field is calculated from. Combining Eqs. ^ 
and 

{Tr{Mp}) = ^a{l + cos0 {cos(j))) + ^(3 {1 ~ cos9 (cos0)), (9) 

where brackets indicate expectation value and average over the different configurations of the nuclear system. In Eq. 
([9]) we assume that the nuclear magnetization is negligible and that _BiV(and therefore 0) has a symmetric distribution 
(Te., l^cfp'^^^^ = 0, A;=l,2,3. . . ). By comparison with the case in which no nuclear field is present and in the limit 
((/)^) < 1, we define the signal Sa as 

Sa ^ {Tr{Mp})^ - {Tr{Mp})^^, 

= (a-/3)cos(?i^^^^^|^e-(^*/^-)"«^(^~a)cos0e-(^*/^-)" ■ ^ ^ 

In deriving Eq. ( |10| wc introduced the coherence decay of the sensor spin characterized by the relaxation time T2e and 
the exponent ptIM] 7 ~ 3. Note that the presence of the nuclear field translates into a change of the NV center average 
fluorescence proportional to the nuclear-spin-induced phase variance. The 'signal' amplitude also grows linearly with 
the difference between the average fluorescence in each of the two possible spin states and reaches a maximum value 
when the phase difference 9 between the excitation and projection pulses is either zero or a multiple of tt (see Fig. 1). 

In order to determine the limiting signal-to-noise ratio, we make use of the property = a*'|0)(0| -I- &'^|1)(1| and 
that Aa^(A&^) = a(/3) for Poisson variables, to calculate the variance 

AAf2 ^ (Tr [M'^p]) - {Tr{Mp}f ^ 

= \{a + l3) + \{a-p)cose (cos 0) e-(^*/^-)" + i (" " Pf {} - cos' ^ (cos e-(A*/T2.)-'^ . (H) 
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The signal-to-noise ratio, SN R=Sa/ is then 



SNR 



-1 _ l\M' 



gg2(At/T2,)^ r (a+;3) + (a-/3) cos 6 e-(At/T2,)T ~l 4^2(At/T,, J l-cos''e(l-(0"))e-'^'/^2e)^ I (12) 



(<^2)2 I (Q-^)2cos2e (02)2 ^ cos2 6 

where, for simpHcity, wc have assumed ((/)^) < 1 and cos^ 6* 7^ 0. Note that in the hmit e~^^*/"^2s)^ ^ 1 the first 
(otherwise dominant) term can be cancelled if we choose — tt and assume that jm, = 1) is a 'dark' state (i.e., /3 = 0). 
The latter, however, is not always the case in practice because, as pointed above, we have a = 1.5/? for direct NV 



spin detection. Therefore, we recast ( 12 ) in the approximate form 

where we made use of the fact that in current experimental settings a = 1.5 x 10~^ <C 1. [SI HO] Hence, the optimal 
sensing time becomes a compromise between the increase in SNR due to larger phase change {(p^) and the exponen- 
tial decay due to decoherence. A similar sensitivity limit is obtained from the measurement of the signal fluctuation, 
as explained in Appendix A. 



Starting from ( 13 1, we can obtain a numerical estimate of the total time T necessary for SNR = 10: At a distance 
d ^ 15 nm from the surface, and for a densely protonated sample we use Q to find AB^ ^ 110 nT. For a sensing 
interval At ^ 40/xs <C ~ 1 ms we get (0^) = 0.6, thus requiring NA^ 2 x 10^ repetitions and a total time 
Tp = NaAI = 60 s (note that in the present case tp^ep, tread At, see Fig. 1). This sensitivity limit could be 
improved enormously if single-shot read out was available. Some strategies toward single-shot readout have recently 
been proposed, such as better collection efficiency via coupling of NV center to a nano-photonic wave-guide [11| or 
readout enhanced by a nuclear spin memory. |12l I13j In this last strategy, nearby nuclear spins (such as the nitrogen 
associated with the NV center or a ^^C) are used to store the information regarding the state of the electronic NV 
spin, so that a given measurement can be repeated many times by mapping back the state of nuclear spin onto the 
electronic spin after each readout. With this technique, it is possible to further improve the SNR although at the 
expense of a much longer readout time (approaching several milliseconds). 



IV. MEASUREMENT OF NUCLEAR SPIN TIME CORRELATIONS 



In the previous section, we implicitly assume that the nuclear correlation time T2n is smaller than the single mea- 
surement time (in practice, of order ^ At) since successive measurements must be independent if they are to improve 
the SNR. However, the opposite regime ^ At allows one to extract valuable spectroscopic information on the 
sample system. Intuitively, this is possible because, as nuclei evolve coherently from a random initial state, the corre- 
lation function — and thus the power spectrum — of sample spins can be determined from the statistics of successive, 
time-delayed measurements. |14| Consistent with the assumption that the nuclear system evolves unperturbed by the 
NV center and that it behaves as a classical magnetic field, we define the autocorrelation function 



Km (t) = {Tr {Mp {At + r)} Tr {Mp [At)}) 



(14) 



with p{t + At) denoting the density matrix that evolved under the action of the nuclear field between the times t and 



t + At (thus acquiring the phase (j>{t + At) = J^~^^* IbBn (t') dt'). Note that since the phase acquisition takes a time 



At, me must restrict t in Eq. (14) and thereafter to r > At. Combining Eqs. (14) and ([2]) we find 

Km (r) 



^ (a + /3) V ^ (a - /3)^ {Kc (t) cos^ 9 + if, (r) sin^ 61 + 2 cos 6* (cos <j))) 



(15) 



where Kc{t) = (cos0(T + At) cos (At)) and Ks{t) = {sin cj) {t + At) sin (At)). Using < 1, and choosing 
9 ~ {2k + 1) 7r/2, we recast Eq. (14) in the simpler form 



^Af(r)-^(a + /3)Vi(a-/3)^ 



'{T + At)(b{At)) 



(16) 



Eqs. (16) and (10) can be used to reconstruct the autocorrelation function K^{t) = {<p{T + At) <p{At)) and to 
determine the sample power spectral density of the phase / — here having the role of a stochastic variable describing 
a stationary random process — via the Wiener-Khintchine theorem |14j 



K^ (r) e- 



'dT. 



(17) 
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Note that because of the finite phase acquisition time, Eq. (17) is restricted to a bandwidth defined by the inverse 
of the separation between two successive measurements 1/At (and has a central observation frequency determined 
by n/ At, with n representing the number of 7r-pulses within the contact time At). 



V. IMAGING AND SPECTROSCOPY OF NUCLEAR SPINS IN BIOLOGICAL SYSTEMS 



In this section we consider examples that illustrate some of the potential advantages — and limitations — of using the 
proposed technique to reconstruct an image or a local nuclear spin spectrum. In each of the simulations that follow 
we use a virtual 'sample spin source' that we recreate in the most realistic way possible from results obtained with 
other techniques. For image reconstruction purposes we assume that the sensor — in the form of a cantilever-mounted 
scanning NV center — can be positioned relative to the sample surface with nanoscale precision. 




FIG. 2; (Upper left insert) High-resolution SEM image of fixated E. Coli. Brighter (darker) regions correlate with high (low) 
spin density. (Main images) 'Raster scan' reconstructions of the corresponding 2D spin lattice. The color code gauges the 
average NV center fiuorescence as determined after NA observations during which the spin alignment changes randomly (see 
Eq. (10)). The darker regions between and on the surface of the bacteria are artifacts resulting from artificial shadowing of 
the source SEM image. The tip-NV center distance and raster scan resolution is 30 nm (left) and 15 nm (right). The scale bar 
corresponds to 300 nm. Other parameters as listed in the text. 



We start by considering the SEM image of Escherichia coli shown in the insert to Fig. 2. Specimens of this 
kind are usually fixated in a dry environment to preserve its morphology meaning that the color code in the image 
correlates with the spin density of protons. The two main images in Fig. 2 show the result of our simulations for 
which we considered the sample as a collection of classical, independent magnetic dipoles with a short correlation time 
(see below). A grey scale is used to indicate the average fluorescence of the NV center at each position [Sa in Eq. 



10). In our calculations, spins were distributed on a regular lattice with 1 nm separation and were given amplitude 



proportional to the local proton density (as implied by the SEM source image) . The NV center distance to the sample 
surface was kept at d =15 nm in one case (right) and 30 nm in the other (left). The evolution time At was 40 and 
120/iS respectively and the number of measurements per image point was Na = 4 x 10^. The resulting time per pixel 
Tp is 15 s (or 50 s) and the image time is estimated at 5 hs (or 4 hs) for a square of (500 nm)^. We note that 
longer exposure times will be necessary if other, non-fundamental sources of noise are present; this scenario, however, 
is unlikely in an optimized confocal microscope where operation has been shown to be photon shot-noise limited. ^ 
One aspect of our example that deserves special consideration concerns the values assumed for the nuclear correlation 
and electron coherence times. First, we note that after fixation the system of Fig. 2 can be considered a solid with 
the result that the nuclear correlation time Tjv -assuming an external field Ba stronger than the internuclear dipolar 
interaction- is dictated by Jiat, the nuclear spin-lattice relaxation time. Because under realistic conditions Ti^r 
largely exceeds At, the time required for a raster scan of the sample grows to impractical values if the nuclear spin 
configuration at a given position must change randomly before the next measurement is carried out. Fortunately, 
there are ways to circumvent this problem, the simplest being to probe other points of the sample surface during the 
wait time. 

Attaining the longest coherence time in a NV center -exceeding 1 ms in isotopically depleted samples [1]- demands 
intercalating a 7r-pulse at the midpoint of the evolution time At. [9, 10 While, for simplicity, our calculations have 
obviated this need [20j , one immediate practical consequence is that a synchronous 7r-rotation must be applied on the 
sample spins if the net effect of the sample dipolar field on the sensor is to be preserved. In the presence of an auxiliary 
dc field Ba ~ 200 gauss, the latter can be carried out via a resonant 'radio-frequency' pulse (at ~ 1 MHzi.pi] 

Although spin density mapping is the most basic form of imaging, it is ultimately the ability to introduce contrast 
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between different soft tissues tliat separates MRI from otlier imaging technologies. Some of the contrast methodologies 
used in MRI find a natural extension in our detection strategy. For example, molecular diffusion away from the 
immediate vicinity of the sensor results in a shortening of the nuclear correlation time, which, with a suitable selection 
of At, can be exploited to make the time-averaged phase shift (p negligibly small. The latter is shown in Fig. 3 where 
we used an SEM image from the membrane of a red blood cell to encode the correlation time of spins on a uniformly 
dense lattice (i.e., spins in the void spaces of the SEM image were assigned a shorter nuclear correlation time). This 
example provides a rudimentary model for a 'water-filled' membrane whose semi-rigid skeleton can be distinguished 
from the embedded fluid. 




FIG. 3: (Insert) SEM image of the membrane of a red blood cell. Void spaces become apparent only after dehydration and 
fixation. (Main images) Unlike Fig. |2] the virtual 2D spin matrix is uniform (emulating the case of a 'wet membrane'). This 
time the color scale of the source image was used to encode the local nuclear spin correlation time. In the example presented 
on the left, 'mobile' regions (corresponding to dark regions in the source image) have a correlation time only 1.3 times shorter 
than the rest. The image on the right is based upon identical conditions except that the correlation time difference was three 
times greater. The scale bar corresponds to 300 nm. 



In situations similar to that of Fig. 3, Eq. (17 1 could be used, for example, to monitor diffusion processes. In this 
context, we note that one of the most important structural characteristics of the cell membrane is that it behaves 
like a two-dimensional liquid, i.e., its constituent molecules rapidly move about in the membrane plane. Therefore, 
one could imagine extensions of the basic pulse protocol to emulate their corresponding NMR counterparts (but with 
resolution on the tens of nanometers). In principle, a broad range of diffusion rates is within reach (because the 
probing time can be greatly enhanced if, after a given evolution period, the NV center coherence is stored in an 
adjacent ^"^C nucleus for future retrieval [9]). Studies of this kind may prove worthy, especially if we keep in mind 
that although the structure of plasma membranes is known to be inhomogeneous, the precise architecture of this 
important system still remains unclear. |15| 

In a different implementation where the auxiliary field Ba points along an axis non-coUinear with the crystal field 
one could rely on the above formalism to extract spectroscopic information from random nuclear spin coherences. An 
example is shown in Fig. 4 where we consider a set of (model) molecules with a 13 Hz heteronuclear (e.g., proton- 
phosphorous) J— coupling. In our simulation the auxiliary magnetic field Ba is 20 gauss, the tip distance is 30 nm, 
and the system correlation time is 100 ms. Assuming the sensor at a fixed position in space. Fig. 4 shows the pulse 
sequence and resulting autocorrelation function and power spectral density. Implicit in this model is the idea that 
molecules tumble and move relative to each other so as to cancel inter- and intramolecular dipolar couplings without 
escaping the observation volume of the sensor during At. Our example mimics the conditions of 'restricted diffusion' 
found, for example, within a cell membrane where molecules 'hop' between adjacent, ~ (lOOnm)"^ compartments on 
a time scale of several milliseconds. [IB] 



VI. CONCLUSION 



While high-field MRI serves as a superb tool to probe the living world, achieving submicroscopic spatial resolution 
presently appears to be a goal exceedingly difficult. Indirect detection via NV centers in diamond provides an alterna- 
tive platform that we examined by means of analytical and numerical calculations. We considered the particular case 
of a single NV center interacting with a large number of nuclear spins, a condition that we described in semiclassical 
terms. When brought in close proximity to the sample surface, e.g. with the aid of a high-precision scanner, the NV 
center is selectively sensitive to field fluctuations induced by nuclear spins immediately adjacent to the sensor (even if 
the mean sample magnetization is negligible). The important practical consequence is that pre-polarization magnet. 
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FIG. 4: In this example the NV center repeatedly monitors a set of equivalent protons subject to a 6.5 Hz heteronuclear 
J-coupling with a second (invisible) spin-1/2 species. Depending on the alignment of the latter, protons precess with one of two 
possible frequencies. (Top) Schematics of the pulse sequence; n denotes the number of p-pulses within the evolution interval 
At. (Bottom) Reconstructed correlation (insert) and corresponding spectral density. Note the factor 2 in the observed splitting 
(13 Hz), a direct consequence of having assumed |cos0| = 1 ('quadratic response'). In the simulation n=l, d =30 nm and 
Ai=100 ms. The nuclear correlation time is 100 ms, and the number of single measurement pairs per point in the correlation 
curve is 4x105. The external magnetic field BaIs 5 mT and pulses acting on nuclear spins are assumed to be broadband so as 
to invert proton spins as well as the J-coupled species. Other conditions as listed in the text. 



gradient coils, and fast-switching current amplifiers — today mandatory in a nuclear spin imaging experiment — are not 
requisites of this technology. 

Our calculations show that simple Ramsey or spin-echo sequences are able to probe the nuclear spin system although 
the relative phase between pulses plays a crucial role. Under current experimental conditions, photon shot noise is the 
main source of error. We stress, however, that the sources of this limitation are not fundamental and that technical 
advances could lead to significant decrease in the imaging times. 

When compared to other kinds of microscopies, several distinguishing features of NV center-based magnetometry 
emerge. For example, given the sharp dependence on the sensor-sample distance, detection is restricted to surface 
spins (with the result that careful sample preparation will be necessary when inner structures of a system are to be 
exposed). On the other hand, the same setup could be exploited to reconstruct 3D topographic maps that can then be 
used to enrich the information content of the images produced via the NV center fluorescence. Even if exposure times 
longer than those typical of other imaging schemes are necessary, diamond-based magnetometry has the potential 
to gauge changes in the dynamics and chemical composition of the sample, thus opening the door to various types 
of contrast. In particular, we have shown that, with an adequate protocol, one could probe molecular diffusion or 
reconstruct the low- or zero- field nuclear spin spectrum [171118] with nanoscale spatial resolution. Finally, we note that 
most biochemical reactions are thermally-driven, stochastic processes that involve the crossing of a barrier or diffusion 
over some kind of potential energy surface. Therefore, the ability to conduct experiments in an open environment, at 
room temperatures can prove crucial to expose the dynamics of living systems in ways not possible with traditional 
magnetic resonance. For example, with spatial resolution of ~ 5 nm — only slightly better than our target here — one 
could envision investigating the stepping of single molecular motors, a process that usually takes place in the tens of 
milliseconds range. 
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Appendix A 



While through Eq. ( 10 1 we monitor sample spins via changes in the average number of photons emitted by the sensor, 
similar information can be obtained if we measure instead changes in the fluorescence variance. This strategy closely 
mimics that already demonstrated in magnetic force microscopy? and, within the framework presented above, appears 



as a 'natural' alternate pathway. Starting from ( 11 ) and neglecting for simplicity relaxation, a simple calculation shows 
that the signal Sy in this case is given by 

S-y = AAf2 - AAf 9^ lcos6' {a - l3) {{a - [3) cos 6 - 1) {(t)^) . (Al) 

To determine the limit uncertainty, we define the auxiliary operator V = {M — {Tr {M p\))^ and calculate /SV"^ — 
(tt ^{V ~ {Tr {V p}))^ . In the limit in which the shot noise is stronger than the spin noise, we find after a 
lengthy but straightforward calculation AV^^ = a. Therefore, the signal to noise ratio is given in this case by 



SNR^^^ 0.1 ^/S, (A2) 



in agreement with ( 13 ) 



We note that Sa in Eq. (10 1 -and thus Sy, Eq. (All- is insensitive to fluctuations of the nuclear field if the 



phase difference 9 between the excitation and projection pulses is an odd multiple of 7r/2. In a way, this condition is 
counterintuitive because, in a sequence where the pulses are phase-shifted, the magnetometer responds linearly — not 
quadratically — to external fluctuations prompting the question as to whether higher sensitivity can be reached. 

Though in a different context, Wincland and collaborators have discussed similar problems extensively. [191 Their 
work highlights the ambiguity that stems from the quantum character of the sensor via the concept of 'quantum 
projection noise': When a single two-level system probes a (non-fiuctuating) magnetic field, maximum sensitivity 
comes at the price of complete uncertainty in the outcome of a measurement; reciprocally, when the measurement 
variance is zero, so is the sensitivity to external fields. Although in the present case the 'signal' comes in the form 
of fluctuations of the magnetometer phase, this principle does play here an important (if more subtle) role. We can 



make it explicit by rewriting (11 1 as 

AAf2 = AAf2 + AAf^, (A3) 

where AAf2 = (tt {iVP p] - {Tr{Mp}f'^ and AAf^ ^ (^{Tr {M p}f'^-{Tr {M p}f . The first contribution measures 

the 'quantum projection noise' or uncertainty in a population measurement of a single two level system; the second 
term corresponds to the nuclear-spin-noise-induced variance in a 'classical', macroscopic- like sensor (where the average 
polarization can be determined from a single measurement). If, for simplicity, we consider in Eq. ([2]) a = 1 and 6 = 0, 
we find after a simple calculation 

AA//2 = i(l-(cos2(0 + 0))), (A4) 



and 



AMI = \ ((cos" + ^)) - (cos (0 + (?))') . (A5) 

For the special case d — {2k + 1)^/2 it follows AM^ — ^ (sin^ 0) and AAI^ = ^ (l — (sin^ 0)) meaning that as we 
increase the amplitude of the external field fluctuations, the gain in the 'classical' contribution to the variance is lost 
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because of an equal but opposite change of the quantum projection noise. This is no longer the case when 9 = kn 
thus leading to an observable change in the variance of the sensor fluorescence. 
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